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ABSTRACT 

We present the detection of the SiO v=Q J=2-l and J=3-2 lines, and of the HCO l O i-0 O o .7=3/2-1/2 
F—2-1 line at ambient velocities towards the molecular outflows in L1448-mm and L1448-IRS3. This 
is the first detection of HCO in a dark cloud. We have also measured lines of H 13 CO + , H 13 CN, HN 13 C, 
CH3OH, and N2H+. While the HCO and the SiO lines have the narrowest profiles with linewidths 
of ~0.5kms _1 , the other lines have widths of ~lkms _1 . Towards L1448-mm, all lines except those 
of SiO and HCO, show two distinct velocity components centered at 4.7 and 5.2 kins" 1 . HCO is 
only observed in the 4.7 km s -1 cloud, and SiO in the 5.2 km s -1 component. The SiO abundance is 
~10 -11 in the 5.2 kms" 1 clouds, one order of magnitude larger than in the 4.7 km s -1 component and 
in other dark clouds. The HCO abundance is ~10~ 11 , similar to that predicted by the ion-molecule 
reactions models for the quiescent gas in dark clouds. The large change in the SiO/HCO abundance 
ratio (>150) from the 4.7 to the 5.2 km s -1 component, and the distribution and kinematics of the SiO 
emission towards L1448-mm suggest that the ambient SiO is associated with the molecular outflows. 
We propose that the narrow linewidths and the abundances of SiO in the ambient gas are produced by 
the interaction of the magnetic and/or radiative precursors of the shocks with the clumpy pre-shocked 
ambient gas. 

Subject headings: stars: formation — ISM: individual (L1448) — ISM: jets and outflows — ISM: 
molecules — ISM: structure 



1. INTRODUCTION 

It is well-known that the SiO emission traces the mate- 
rial which has been processed by shocks. The SiO profiles 
present very broad linewidths that are centered at veloc- 
ities different from that of the ambient cloud (Martin- 
Pintado, Bachiller, & Fuente 1992). As silicon is heav- 
ily depleted onto dust grains in the cold quiescent gas, 
the SiO abundances are extremely low in dark clouds 
(<1CT 12 ; Ziurys, Friberg, & Irvine 1989; Marti'n-Pintado 
et al. 1992). Nevertheless, Lefloch et al. (1998) reported 
the first detection of bright and narrow (with the same 
linewidth of the ambient gas) SiO lines at ambient veloc- 
ities associated with the molecular outflows in the dark 
cloud NGC 1333 (t>10 4 years). Lefloch et al. (1998) sug- 
gested that the SiO emission arises from the postshock 
equilibrium gas after the interaction of powerful proto- 
stellar jets with local dense clumps. They concluded that 
the narrow SiO emission could be only detected in partic- 
ular objects characterized by its high degree of dumpi- 
ness. Codella, Bachiller, & Reipurth (1999) confirmed 
the suggestion of Lefloch et al. (1998) that the ambient 
SiO emission arises from decelerated postshock material. 
In this scenario one would not expect to detect SiO ambi- 
ent emission with narrow linewidths in very young molec- 
ular outflows. 
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The L1448 outflow is a very young outflow with a dy- 
namical age of only ~3500 years (Bachiller et al. 1990). 
Due its youth, this object provides an excellent labora- 
tory to study the origin of the ambient SiO emission, 
because the time scale is too short for the postshock gas 
to reach the equilibrium (Chieze, Pineau des Forets, & 
Flower 1998) and to decelerate to the radial velocities of 
those of the ambient gas. 

The formyl radical (HCO) has not been detected so far 
in dark clouds (e.g. L183, where the derived HCO abun- 
dance is <2xl0 -10 ; Schenewerk et al. 1988). However, 
gas-phase chemistry models predict even lower HCO 
abundances (^KU 11 ) than the previously measured up- 
per limits to this molecule (Leung, Herbst, & Hueb- 
ner 1984). 

In this Letter we present the first detection of very nar- 
row SiO and HCO emissions at ambient velocities in the 
L1448-mm/IRS3 molecular outflows. The narrow lines 
and the radial velocities of the SiO emission do not seem 
to be related to decelerated postshock equilibrium gas 
but produced by the interaction of the magnetic and/or 
radiative shock-precursors with the ambient pre-shocked 
clumpy medium. 

2. OBSERVATIONS AND RESULTS 

The observations were carried out with the IRAM 
30-m telescope at Pico Veleta (Spain). We observed 
the central position of L1448-mm (the driving source 
of the highly collimated outflow) and L1448-IRS3, and 
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Fig. 1. — Sample of line profiles for the different positions to- 
wards the L1448-mm/IRS3 molecular outflows. The offsets, in 
arcseconds at the top of the panels, refer to the position of L1448- 
mm (a(2000) = 03 /l 25 m 38*0, 5(2000) = 30°44'05"). The dashed 
vertical lines show the radial velocities of the velocity components 
for which SiO and HCO have been detected. 



mapped the L1448-mm region observing towards the 
offsets (20,0), (-20,0), (0,20), and (0,-20). Table 1 
summarizes the parameters of the observed transitions. 
The beam size of the telescope was ~27" at ~90 GHz 
and -18" at ~140 GHz. The 3 mm and the 2 mm 
SIS receivers were tuned to single side band with im- 
age rejections >10 dB. We observed simultaneously the 
SiO J=2-l, H 13 CO+, HN 13 C, and HCO species using 
the frequency-switching mode with a frequency throw of 
3 MHz. As spectrometers we used VESPA, which pro- 
vided a velocity resolution of ~0.14kms . Typical sys- 
tem temperatures were 90-100 K. We observed the rest of 
the species (H 13 CN, CH 3 OH, N 2 H+, and the SiO J=3- 
2 transition) in position-switching mode. The system 
temperatures and the spectral resolutions were typically 
of 120-350 K and 0.13-0.20 km s" 1 , respectively. All the 
line intensities were calibrated in antenna temperature. 

Fig. 1 shows all the spectra measured for the L1448- 
mm map and L1448-IRS3, and Table 1 shows the ob- 
served parameters for all the lines measured towards 
the central positions in L1448-mm and L1448-IRS3. 
For comparison, Fig. 1 shows, in the upper panels, 
the profiles of the blueshifted and redshifted shocked 
gas with broad, — 50kms _1 , lines (see e.g., Martin- 
Pintado et al. 1992). The ambient profiles are shown 
in detail in the central panels of Fig. 1. The SiO and 
HCO ambient emissions show very narrow linewidths of 
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Fig. 2. — Correlation between the narrow SiO emission inten- 
sity and the ion (H 13 CO + ) to neutral (HN 13 C) intensity ratio for 
the 5.2 and 4.7kms velocity components detected towards the 
positions mapped in L1448-mm (the position offsets shown above 
the points are relative to L1448— mm). The vertical arrows indicate 
upper limits to the SiO emission. Note the trend of the relative ion 
to neutral abundance to increase at the positions where SiO has 
been detected. 



— 0.5kms , smaller than those of the ambient cloud 
molecules like H 13 CO+, HN 13 C, and N 2 H+ which have 
widths of ~lkms _1 . There are appreciable differences 
between the observed profiles of SiO, HCO, and the other 
molecules. While all ambient molecules towards L1448- 
mm show double peaked profiles with velocity peaks at 
4.7 and 5.2 km s -1 (see vertical dashed lines in Fig. 1), 
the HCO and SiO lines are single gaussian narrow pro- 
files with peaks at 4.7 and 5.2 km s" 1 respectively (Ta- 
ble 1). VLA images of the ambient NH3 gas show that 
these two components correspond to distinct molecular 
clumps with different spatial distribution associated with 
the L1448-mm source (Curiel et al. 1999). 

From our limited map, we notice that the detection 
of the SiO emission in the 5.2 km s -1 cloud is related to 
an increase of the ion abundance with respect to neu- 
trals. Fig. 2 illustrates this effect by representing the 
SiO line intensities as a function of the ratio between 
the ion (H 13 CO+) and the neutral (HN 13 C) intensities 
for the 4.7 and the 5.2 kms -1 components. Although 
the data are scarce, the tendency of the SiO emission to 
increase with the ion to neutral intensity ratio is clear. 
Since the lines are optically thin, this shows that the 
larger SiO abundance in the ambient cloud is associated 
with an enhancement of the ion abundance relative to 
the neutrals. Furthermore, we find that the SiO profile 
towards L1448-mm (0,-20) presents a substantial broad- 
ening (redshifted with respect to the ambient gas) due 
to material which has already entered the shock. Al- 
though CH3OH is also a neutral molecule, its emission 
is shifted to the 5.2 km s -1 component. As discussed in 
section 4, these data are consistent with the fact that 
the SiO and CH3OH abundances are enhanced by grain 
mantle destruction produced by shocks (Millar, Herbst, 
& Charnley 1991). 

Towards L1448-IRS3, we have measured the strongest 
narrow SiO emission in L1448. Again, the SiO and HCO 
lines are narrower (with linewidths — 0.5kms _1 ) than 
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a Rcfcrs to the F = — + 1 hypcrfinc component. 
b Refers to the J = 3/2 — > 1/2, F = 2 -> 1 transition. 
c Rcfcrs to the F = 2 — » 1 hypcrfinc component. 

d Intcnsities are in antenna temperature. The limits correspond to 3cr. 



those of the ambient molecules. Toward this position 
all line profiles are single peaked, but the SiO emission 
shows an appreciable velocity shift relative to the emis- 
sions from ambient molecules (~0.2kms _1 ; see vertical 
line in Fig. 1). Like in L1448-mm, this narrow emission 
is shifted in the same velocity direction as the shocked 
material which suggests that the "ambient" SiO is re- 
lated to shocks. 

3. ABUNDANCES 

Assuming a kinetic temperature of 20 K (Curiel et 
al. 1999), Table 2 shows the derived H2 densities and 
the column densities of all molecules for the central re- 
gion of L1448-mm and L1448-IRS3. The H 2 densities 
derived from the SiO J— 2-1 and J=3 -2 lines are a few 
10 5 -10 6 cm~ 3 . If we consider a HCO + fractional abun- 
dance of ~1CT 8 (Irvine et al. 1987), the derived SiO 
abundance towards L1448-mm (0,0), (-20,0), and (0,- 
20) is ~10~ n for the 5.2 kms^ 1 component, an order of 
magnitude lower than what Lefloch et al. (1998) found in 
NGC 1333. The estimated upper limit to the SiO abun- 
dance for the 4.7kms _1 cloud is <10~ 12 at L1448-mm 
(0,0), of the same order as that obtained in the cold qui- 
escent clouds TMC1 and Bl (Ziurys et al. 1989; Martin- 
Pintado et al. 1992). Therefore, the SiO abundance in 
the 5.2 km s -1 component has been enhanced by more 
than one order of magnitude with respect to the quies- 
cent gas of the 4.7 kms -1 cloud. 

On the other hand, HCO has not been detected for the 
5.2kms _1 component (Fig. 1), with an upper limit of 
x(HCO)<4xl0~ 12 (Table 2). For the 4.7km S - 1 cloud, 
we derive a HCO abundance of few 10 -11 . 

4. DISCUSSION 

HCO is predicted to be produced in dark clouds 
through ion-molecule reactions with abundances of 
lO" 11 (Leung et al. 1984). Our detection of HCO and 
the abundances derived for this molecule (~10 -11 ; Ta- 
ble 2) are in agreement with these models, suggesting 
that the 4.7 km s^ 1 material traces the quiescent gas. 



The SiO/HCO abundance ratio changes by more than 
two orders of magnitude from the 4.7 to the 5.2 kms -1 
component of the ambient gas, indicating that the for- 
mation of both molecules are related to different mecha- 
nisms. 

The spatial distribution and the kinematics of the nar- 
row SiO emission suggest that this emission is associ- 
ated with the molecular outflows towards L1448-mm and 
L1448-IRS3. However, the lines are very narrow and 
with the radial velocities of the ambient gas. We can 
rule out the idea that the ambient SiO emission is gen- 
erated in the postshock equilibrium gas since the time 
scales for the shocked gas to decelerate to the ambient 
velocities are much larger (>10 4 years) than the dynami- 
cal age, and the SiO has not been detected in the ambient 
4.7 kms -1 cloud. 

In the following, we explore the possibility that the 
SiO ambient emission is a signature of the shock precur- 
sors associated with the jets driven by the young stars. 
There is a tendency for the ambient SiO emission to arise 
only from regions where the relative abundance of the 
ions to neutrals is larger than in the quiescent clumps 
(see Fig. 1). This indicates that the ions could have 
slipped to the 5.2 kms -1 component from the neutrals 
of the 4.7 kms -1 cloud, as expected if this narrow emis- 
sion were produced by the magnetic precursor associated 
with C-type shocks. The magnetic precursor would force 
the plasma to stream through the neutral gas ahead the 
jump front, generating a velocity decoupling between the 
charged and neutral fluids (Draine 1980). As the charged 
particles were compressed in the first stages of the mag- 
netic precursor, the ion density would initially increase in 
the 5.2 km s -1 component with respect to the quiescent 
4.7 kms -1 cloud (see e.g. Flower et al. 1996). 

Pilipp, Hartquist, & Havnes (1990) and Pilipp & 
Hartquist (1994) have considered grains as a separated 
fluid from the ion, neutral, and electron components. 
The behavior of the grains changes depending on several 
parameters like the H2 density, mean charge, and grain 
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TABLE 2 
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0.3 


0.4 


2.8 


1.7 




40 


lxlO -11 


2X10" 11 



Note. — Abundance upper limits were estimated from the 3er level. These parameters have been derived for the central region in each source. 



size distribution. There are cases where the charged 
grains are coupled to the ions, leading to a small grain- 
neutral velocity drift just before the shock starts. The 
Chieze et al. (1998) time-dependent C-shock model 
shows that, even for time scales of <10 4 years, this ve- 
locity drift generates the sufficient friction to produce 
the ejection of small amounts of silicon from grains into 
the gas phase in the 5.2 kms -1 cloud through sputtering 
(Flower et al. 1996; Caselli, Hartquist, & Havnes 1997). 
Since these time scales are rather short, we should ex- 
pect that silicon is released directly in the form of SiO, 
producing an abundance of 10~ 10 or 10~ n (Flower et 
al. 1996; Caselli et al. 1997). The lack of a velocity shift 
between the SiO and the ions for the 5.2 km s" 1 compo- 
nent also indicates that we are observing the very early 
phases of the interaction of the shock precursor. This is 
supported by the data at the position (0,-20) where the 
broader SiO profile (see Fig. 1) suggests that in this case 
a large fraction of the neutrals has already entered the 
shock. 

On the other hand, the detection of [Si II] 34.8 [im 
emission in L1448-IRS3 (Nisini et al. 2000) shows the 
existence of a photon-dominated region (PDR) which 
could be related to the radiative precursor of the J- 
type shocks (Shull & McKee 1979; Hollenbach & Mc- 
Kee 1989). In fact, the C-type shock models of Chieze 
et al. (1998) predict a jump discontinuity in the neutral 
flow in the early stages for time scales <10 4 years. The 
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UV radiation field of the radiative precursor associated 
with such J-shock could affect the chemistry of the pre- 
shocked medium. It might be possible that the ambient 
SiO emission would be produced by the radiative pre- 
cursor through the photodesorption of Si from the grains 
(Schilke et al. 2001). However, the SiO/HCO abundance 
ratio in the 5.2kms _1 velocity component is very differ- 
ent from what is observed in PDRs (Schilke et al. 2001). 

In summary, the linewidths, velocities, abundances, 
and spatial distribution of the SiO in L1448-mm/IRS3 
could be explained by the radiative and/or magnetic 
shock precursors associated with the jets driven by pro- 
tostars. Based on all these results, we speculate that we 
have probably detected the first evidence of the chemical 
effects on the ambient cloud produced by the precursors 
of interstellar shocks. Higher angular resolution obser- 
vations of SiO, HCO, and tracers of the ambient gas are 
needed to clearly establish the origin of the ambient emis- 
sion of SiO in L1448. 
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